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We study a two-dimensional motion of a charged particle in a weak random potential and a 
perpendicular magnetic field. The correlation length of the potential is assumed to be much larger 
than the de Broglie wavelength. Under such conditions, the motion on not too large length scales 
is described by classical equations of motion. We show that the phase-space averaged diffusion 
coefficient is given by Drude-Lorentz formula only at magnetic fields B smaller than certain value Be. 
At larger fields, the chaotic motion is suppressed and the diffusion coefficient becomes exponentially 
small. In addition, we calculate the quantum-mechanical localization length as a function of B in 
the minima of axx- At B < Be it is exponentially large but decreases with increasing B. At B > Be, 
the localization length drops precipitously, and ceases to be exponentially large at a field B* , which 
is only slightly above Be- Implications for the crossover from the Shubnikov-de Haas oscillations to 
the quantum Hall effect are discussed. 



I. INTRODUCTION 

In this paper we study a two-dimensional motion of a 
charged particle in a weak random potential and a per- 
pendicular magnetic field. This problem has deep his- 
torical roots and the limiting cases of a weak and a very 
strong magnetic field are fairly well understood. As we 
will see below, the nature of the motion in these two lim- 
its is crucially different. Surprisingly, until now no theory 
for the cross-over between the two limits has been pro- 
posed. Our goal is to develop such a theory. We will 
start with a classical description of the transport. 

An important prediction of the classical magnetotrans- 
port theory is that the conductivity in the direction per- 
pendicular to the magnetic field is reduced. 
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where ctq is the zero field conductivity (the magnetic field 
B is assumed to be along the z-direction), ujc — eB/mc 
is the cyclotron frequency, and r is the transport time 
determined by the properties of the random potential. 
Strictly speaking, in classical theory it is more consistent 
to study the diffusion c oeffi cient D. So, we would write 
Drude-Lorentz formula (1.1) in the form 
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where Dq — ^v t is the diffusion coefficient in zero field, 
V being the particle velocity. Drude-Lorentz formula pre- 
dicts that if the magnetic field is not too weak so that 
u!cT > 1, then the diffusion coefficient falls off inversely 
proportional to the square of the magnetic field. 

Let us examine the physical picture of the motion in 
such magnetic fields. It is easy to verify that the Lorentz 



force has a dominant effect on the motion and the devia- 
tions from the perfectly circular cyclotron orbit are small. 
In such circumstances, the original coordinates r — {x, y) 
are not very useful any more. Instead, it is convenient 
to study the motion of the guiding center p = (p^, Py) of 
the cyclotron orbit. 

Suppose the cyclotron gyration is clockwise (this is the 
case if, e.g., the particle charge is negative and the mag- 
netic field is in the negative i:-direction) . The guiding 
center coordinates are defined as follows, 
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Drude-Lorentz formula (1.2) results from the assumption 
that the guiding center p performs a random walk. The 
characteristic step of such a random walk is the cyclotron 
radius. Re — v/lUc, and the time interval between the 
steps is the transport time r. As we will see below, this 
is the correct description of the motion if the magnetic 
field is not too strong. 

Perhaps the first work that demonstrated that Drude- 
Lorentz formula may not be valid, in the limit of strong 
magnetic field was that of Alfventl where he studied the 
motion of a charged particle in an inhomogeneous elec- 
tromagnetic field. This and subsequent studya Q has led 
to the recognition that instead of the random walk, the 
guiding center performs a slow adiabatic drift along some 
well defined contours. The attention to this problem was 
stimulated by its plasma physics applications, and mostly 
three-dimensional case was considered. Not so long ago, 
the extension to the twctdimensional case has been pro- 
posed by several authorsQ motivated by the quantum Hall 
effect studies.Q We will discuss the two-dimensional case 
from now on. 

Conventionally, the drift approximation is applied to 
the regime where the magnetic fields so strong that the 



cyclotron radius Re = v/ujc is smaller than the correla- 
tion length d of the random potential. In this case the 
guiding center performs a drift along the constant energy 
contours of the random potential. For the potential of a 
general type all of such contours except one are closed 
loops and thus the motion is finite. The motion is infi- 
nite only when a guiding center happens to be on the so- 
called percolating contour Jj If one takes the drift picture 
literally, and attempts to calculate the average diffusion 
coefficient, the result will be equal to zero because the 
percolating contour has zero measure. 

Certainly, it has been understood that the drift pic- 
ture is only an approximation. Nevertheless, the diffu- 
sion coefficient should be significantly smaller than the 
Drude-Lorentz result (|l.2| ). We will show that the diffu- 
sion coefficient is, in fact, exponentially small. 

Comparing the transport properties in the two regimes 
described above, we see that the increase in the magnetic 
field drives the system from the essentially delocalized, 
chaotic regime to the regime where the motion is regular 
and the trajectories of the particles are localized. We call 
this phenomenon the "classical localization." The clas- 
sical localization occurs because of an extremely ineffec- 
tive energy exchange between two degrees of freedom of 
the particle, the cyclotron motion and the guiding center 
motion. Without such an exchange the guiding center 
is bound to a certain constant energy contour. At the 
same time, the energy exchange is suppressed because 
the two degrees of freedom have very different character- 
istic frequencies, the cyclotron frequency tOc being much 
larger than the drift frequency cod- Naturally, the present 
problem is directly related to the problem of a nonconser- 
vation of adiabatic invariants. The latter is known to be 
exponentially smallJZI and therefore, it is not so surprising 
that the diffusion coefficient turns out to be exponentially 
small as well. 

One of the quantities we calculate in this paper is the 
value Be of the magnetic field where the diffusion gives in 
to the classical localization as B increases. A naive guess 
would be the field where Re = d. Let us, however, com- 
pare ujc and ujd at such a field. Denote the amplitude 
of the random potential U(r) by W. We will assume 
that the potential is weak, W <^ E, where E = mv'^/2 
is the particle's energy. The characteristic drift velocity 
is Vd ^ VU /mujc ~ W/muJcd, and the drift frequency 
Wd ^ Vd/d '^ W/muJcd^. Hence, the ratio 7 of the two 
frequencies is 
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We see that at the point where Re — d, this ratio is of 
the order of W/E ^ 1. Surprisingly, the classical local- 
ization must first arise already when Re ^ d. To under- 
stand what kind of drift takes place in this case one can 
use the averaging method. This method was extensively 
developed by Krylov, Bogolyubov, and Mitropolskya and 
in application to the problem at hand by Kruskal.H In the 



spirit of this method, one has to imagine that the slowly 
moving guiding center is entirely "frozen" on the time 
scale of the cyclotron period. One then calculates the 
average potential 
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acting on the particle during one cyclotron rotation. Ac- 
cording to the averaging method, the drift of the guiding 
center is performed along the constant energy contours of 
the averaged potential Uq{p). Jhis conclusion was pre- 
viously reached by Laikhtman.0 If Re ^ d, the average 
potential coincides with the bare one and so, in agree- 
ment with the previous studies, the drift is performed 
along the constant energy contours of the bare potential. 
However, if Re ^ d, then Uq differs from U . The aver- 
aging reduces the amplitude of the potential by a factor 
\J Rejd^ which is the square root of the number of uncor- 
related "cells" of size d along the cyclotron orbit of len gth 
27ri?c- Hence, C/q has the amplitude Wq ^ W \J djRe. 

Now we can find the true boundary Be of the classical 
local izat ion. To this end we have to replace W by Wq in 
Eq. (|l.4[), which gives 
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and then solve 7 = 1 for S. The result is 
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The change of the transport resirne at such field wag 
predicted earlier by Baskin et a/.E2l and by Laikhtman.EI 
These authors noted that the displacement 5r of the guid- 
ing center after one cyclotron period is a decreasing func- 
tion of the magnetic field, 5r ~ ^d in our notations. 
Thus, at B > Be where 7 < 1, such a displacement is 
smaller than the correlation length of the random poten- 
tial. As a result, the scattering by the potential is no 
longer a sequence of uncorrelated acts and the motion 
of the guiding center is different from the random walk, 
which invalidates Eq. (1.2). 

Although the cross-over point Be has been identified 
correctly, the understanding of the transport regime at 
larger magnetic fields remained not entirely satisfactory. 
For example, Baskin et al.cB arrived at a strange con- 
clusion that aX B > Be the diffusion coefficient becomes 
larger than that given by Drude-Lorentz formula (B_j). 
On the other hand, the calculation of Laikhtmarffre- 
lies on the existence of the random inelastic scattering 
processes. In this paper we address the question of zero 
temperature transport where all the scattering acts are 
due to the static random potential only. 

The key point of our approach is that the drift picture 
is albeit excellent but an approximation. A more accu- 
rate analysis given in Sec. || reveals that the diffusion 



occurs not only when the guiding center is situated pre- 
cisely on the percolating contour but also-within a strip 
of finite width, so-called stochastic layer ,llil surrounding 
this contour. Such a layer turns out to be exponentially 
narrow if the magnetic field is larger than Be- As a re- 
sult, the phase-space averaged diffusion coefficient D is 
also exponentially small. 
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Thus, the "classical localization" above Be causes strong 
deviations from the conventional Drude-Lorentz for- 



mula (|f.2|). 

The existence of the stochastic layer around the per- 
colating contour is quite natural. Indeed, the classical 
localization is owing to the fact that drift trajectories 
are closed loops. It turns out that the drift along the 
loops passing sufficiently close to the saddle points of 
the random potential is unstable. The instability is re- 
alized as a slow diffusion of the guiding center in the 
direction transverse to the drift velocity. Suppose that 
the percolation level is Uq — 0. By virtue of a small 
transverse displacement, the particle drifting along the 
contour Uq = —e can move to another closed contour 
Uq = +e. Although this displacement may be small, it 
will, in fact, lead to a much larger displacement at a later 
time because the center of the other loop is typically lo- 
cated a large distance away. Eventually, the particle can 
travel infinitely far from its initial position. This is the 
nature of the diffusion mechanism inside the stochastic 
layer. This mechanism is guaranteed to exist because the 
percolating contour necessarily passes through the saddle 
points. 

The suppression of chaotic motion with increasing 
magnetic field proceeds as follows. At B < Be the chaotic 
motion takes place in the majority of the phase space, 
while the the regular motion is restricted to small stabil- 
ity islands .til In this regime the correlations among the 
scattering acts can be ignored and Eq. (f.l) applies. As 
the magnetic field increases, the regions of regular mo- 
tion expand while the stochastic layer shrinks. Above 
Be the width of the stochastic layer starts to decrease 



exponentially leading to formula (1.8) 



So far, we have discussed a purely classical dynamics. 
One can also study the transport properties of a non- 
interacting electron system quantum-mechanically. Due 
to quantum interference, the conductivity p£ such a sys- 
tem turns out to be length-scale-dependent .lli The knowl- 
edge of classical dynamics enables one to find "classical" 
axx, i-e., the conductivity, which would be measured on 
not too large length scales where effects of quantum inter- 
ference are weak. Classical axx is calculated as a product 
of the classical diffusion coefficient D and the quantum 
density of st ates m/Trfi . It is given by Drude-Lorentz 
formula [lA) aX B < Be and by a different formula 
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aX B > Be- Here kp = j^yhnE is the Fermi wavevector 
{E has the meaning of the Fermi energy). The depen- 
dence of classical axx on B is illustrated graphically by 
Fig. 1^. As one can see, classical Uxx quickly drops above 
B = Be- In Fig. H we indicated one special value of the 
magnetic field, _B*, at which classical Uxx reaches ^, 



B^ ^ BelnG. 
Here we assume that G ^ I, i.e., that 

d » fc/ (§ 



(I.IO) 



(l.II) 



As we will see below this value of the magnetic field 
plays an important role in the quantum transport. 
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FIG. 1. Classical conductivity axx (solid line) and the 
localization length ^o at the QHE conductivity minima 
(solid line with dots) as functions of the magnetic field 
(schematically). Dots serve as a reminder that ^o is de- 
fined at discreet values of the magnetic field. The curves 
are labeled by the equation numbers, which render their 
functional form in the corresponding intervals. 

At this point we would like to remind the reader that 
the true axx, i-e., the one which is measured experimen- 
tally, is the conductivity on a large length scale (of the 
order of the sample size) . The calculation of this quan- 
tity is much more difficult. Similar to the classical trans- 
port theory, there exist two mutually contradicting ap- 
proaches. One is the theory of the Shubnikov-de Haas 
(SdH) effect, which aspires to predict the behavior of 
axx in weak magnetic fields. The other is the theory of 
the quantum Hall effect (QHE), which is conventionally 
applied to strong fields. 

At present, the transition from the SdH regime to the 
QHE is not well understood even for a non-interacting 
system. The traditional explanation of the QHE is based 
on the idea of localization; viz., it is believed that at zero 
temperature an electron can propagate diffusively only if 



its energy is precisely at the center of a Landau level (in 
strong fields). Q This leads to isolated peaks in axx^ which 
are the signature of the jQiffi'- On the other hand, in the 
theory of the SdH effect ,E£E3 the suppression of Oxx is re- 
lated merely to the dips in the density of states between 
neighboring Landau levels, while the idea of localization 
is totally discarded. This crucial difference leads to dif- 
ferent predictions for the conductivity minima. Arguing 
from the QHE standpoint, one expects zero dissipative 
conductivity, whereas the theory of SdH effect predicts a 
finite one. 

In this paper we will advocate the following way to 
resolve this apparent contradiction. We will argue that 
at the QHE conductivity minima the states at the Fermi 
level are localized. At B < B^ where B^ is given by 



Eq. (1.10), the localization length ^o of such states is ex- 
ponentially large but decreases from one minima to the 
next as B increases. Above Be the fall-off of ^o is ex- 
tremely sharp and at B ~ B^, which is only logarithmi- 
cally larger than B^, the localization length stops being 
exponentially large. Consequently, B = B^, is the small- 
est magnetic field at which the observability of the QHE 
does not require exponentially small temperatures. This 
fact motivate us to identify the field B = B^, a.s the start- 
ing point of the QHE. In other words, this is the position 
of the "first" QHE plateau. 

To avoid confusion let us further elaborate on this 
issue. Precisely at zero temperature one will observe 
the QHE peaks. Between the peaks axx will be ex- 
actly zero because of the quantum localization. At fi- 
nite temperature T > inelastic processes appear, which 
break the quantum coherence on length scales exceed- 
ing some temperature-dependent length L^{T). Thus, if 
^0 > L^(T)^ then the quantum localization is not im- 
portant and the QHE features disappear. It is believed, 
that the dependence of L^ on T is some power law.EJ 
Therefore, if ^o is exponentially large, then the inequal- 
ity Co > LifyiT) is met already at exponentially small 
temperatures. 

There is yet another way to see why the observability 
of the QHE require small T when ^q is large. It is known 
from experiment (see the bibliography of Ref. Hq) that 
the low-temperature magnetotransport data at the axx 
minima is consistent with the law 
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which can be interpretedt3 in terms of the variaUe-range 
hopping in the presence of the Coulomb gap.llZI In this 
theory Tq is directly related to ^o. 
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where k, is the dielectric constant of the medium. Deep 
minima of axx are observable only if T <C To- Thus, 
if Co is exponentially large, then the QHE can be ob- 
served only at exponentially small T. So, we reiterate 
once more that in practical terms there exists a starting 



point of the QHE. The precipitous drop of Co(-B) above 
Be leaves only a minimal ambiguity in identifying this 
point with B ^ B^. 

Our calculation of the localization length Co at, 
QHE minima of axx is based on the following ansatz\ 
which we discuss in more detail in Sec. IV, 
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Here go — -^axx is the dimensionless classical conduc- 



tance. Using Eqs. (1.1) and ( |l.9| ), we immediately find 
Co (X exp fc^ll) , Be {WIEY < B < B^, (1.15) 
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The low- field end of the interval in Eq. (1.15) corresponds 

to OJrT ^ 1. 



As one can see from Eqs. (1.15) and (1.16), the local- 
ization length indeed drops precipitously above B = Be- 
At i? = _B*, which is only logarithmically larger than 
Be, go becomes of the order of unity and Co ceases to be 
exponentially large. The dependence of Co on B in the 
interval Be (W/E)^ < B < B^ is illustrated by Fig. ||. 
The dependence of Co on B at even stronger magnetic 
fields, B > B^,, will be discussed in a forthcoming pa- 
per. At this point we can only say that at such fields 
the localization length is determined mainly by quantum 
tunneling and exhibits a power law dependence on B. 

In order to verify our predictions concerning Co (5) ex- 
perimentally, one has to measure axx at ve ry low tem- 
peratures and fit the data to the form (1.12). From such 
a fit one can deduce To, which is directly related to Coj 



see Eq. ( |1.13D . 

The paper is organized as follows. In Sec. we dis- 
cuss the classical dynamics in strong (B ^ Be) magnetic 
fields and demonstrate that the diffusion coefficient is ex- 



ponentially small. In Sec. [II we analyze the same prob- 
lem from the quantum-mechanical point of view. Sec. IV 
is devoted to the calculation of the quantum localization 
length both in strong and weak magnetic fields. Finally, 
in Sec. M we summarize our findings and discuss their 
relation to the experiment. 



II. CLASSICAL DYNAMICS IN STRONG 
MAGNETIC FIELDS 

In this Section we study the classical dynamics of the 
system with the Hamiltonian 
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U{r), A={0,-Bx,0). (2.1) 



It corresponds to a particle with negative charge — e and 
the magnetic field in the negative z-direction. Thus, the 
cyclotron gyration is clockwise. By means of the canon- 
ical transformation with the generating function 



F{x,y,9,py) = muJc 
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we obtain new momenta —dF/dpy = mujcPx and 
—dF/d O = I . In terms of the new variables, the Hamil- 
tonian (^.l[) acquires the following form 



H = lojc + U{px + R cos 6,py ~ R sin ( 
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It is easy t o sec that the pair {px, py) matches the earlier 
definition ( |1.3| ) of the guiding center coordinates. The 
geometrical meaning of the other variables is illustrated 
by Fig. 0. 
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FIG. 2. The guiding center and cyclotron motion co- 
ordinates. 

The equations of motion are 

I dU . \ dU ,^^, 

P^ = ^ir~Tr~^ Py = ir~^r^ (2.3) 

(2.4) 

This system contains four dynamical variables, which 
makes its solution difficult. We can eliminate one of the 
variables, e.g., /, using the energy conservation. To this 
end we need to solve the equation 

S = /wc + U{p, 6, 1) 

for /, or equivalently, the equation 
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for R. For the potential U of an arbitrary strength this 
can be quite cumbersome. However, at least when the 
amplitude W of potential U is small enough, 



a 
it is sufficient to use an approximate solution 
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Condition (2.5) guarantees that the deviation of R from 
Re is much smaller than the correlation length d of po- 
tential U . Under this condition we can al so n eglect the 
deviation of 6 from ojc- As a result, Eqs. (2.3) and (2.4) 
can be treated as the equations of motion for the time- 
dependent Hamiltonian 



H — U{px + Re coscjct, Px — Re sinWcO 



(2.6) 



with Py being the canonical coordinate and mtUcPx be- 
ing the canonical momentum. It is customary to classify 
the systems of this kind as systems with 1^ degrees of 
freedom. 

It is useful to expand Hamiltonian ( ^.6[ ) in the Fourier 
series. 
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with the expansion coefficients given by 
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[compare with Eq. (1.5)]. The new equation of motion 
for Px is 
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and similarly for py. If we drop the sum on the right- 



hand side of Eq. (2^), then the remaining term will de- 
scribe the drift of the guiding center along the contours 
of constant C/q- Such a drift leads to the classical local- 
ization described in the previous Section. The character- 
istic drift frequency is of the order of ujd ^ Wo/moj'^(P, 
where Wq is the amplitude of Uq (see Sec. |). Thus, if 
the parameter 7 = ujd/ujc is small, then all the terms in 



the sum on the right-hand side Eq. (2.9) have frequencies 



much larger than the cod- They can be considered a high- 
frequency perturbation imposed on the "unperturbed" 
drift motion. 

The presence of a small parameter calls for the pertur- 
bation theory treatment (averaging method) developed in 
Refs. |2|-Q. Unfortunately, it is not possible to calculate 
the diffusion coefficient perturbatively because the per- 
turbation theory series converge only asymptotically, i.e., 
they formally diverge for any finite 7. The calculation 
of the diffusion coefRcient requires a different approach 
based on the consideration of the chaotic dynamics of the 
system within a narrow stochastic web surrounding the 
percolating contour of the potential Uo{p). 

Due to an extreme difficulty of the problem, we re- 
strict our consideration by two particular examples: a 
chessboard potential and a Gaussian random potential. 



A. Chessboard geometry 



Consider a chessboard potential 
Uix,y) 
In this case Uq is given by 



I X y 

W cos — + cos - 

\ a a 



Py 



Uo = ~WJo (Rc/d) cos ^ + cos ^ 



(2.10) 



More generally, 
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where jTfc's are the Bessel functions. 

As explained above, one can introduce the dimension- 
less param eter 7, which governs the classical dynamics. 
Equation (2.10) suggests that the appropriate definition 
for 7 is 
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Note that with this definition 7 vanishes whenever 
Rc/d coincides with a zero of Jq. This property is a 
peculiarity of the periodic geometry and it leads to os- 
cillations in the diffusion coefficient with the magnetic 
field, which are well inewn to exist both from theory 
and from experiment J13'E2I This behavior is nonuniversal 
and therefore, is not of primary interest to us. In the fol- 
lowing we will assume that the ratio Rc/d is always close 
to midpoints between the successive zeros of JTq. In this 
case, the dependence of 7 on R^ is given by Eqs. ( |l.4| ) 
and (L6). We will focus on the case 7 <C 1. 

The "unperturbed" motion is described by the Hamil- 
tonian 

Ha = Uoip), 

which is time-independent. Hence, Uq is the integral of 
motion in agreement with the statement that the drift 
is performed along the contours Uq = const. The mo- 
tion has a periodic array of hyperbolic (or saddle) points. 
Some of them, (ttc?, 0), (0,7rd), (— Trd, 0), (0, — Trd) are 
shown in Fig. 0, the others can be obtained by periodic 
translations. The hyperbolic points are connected by het- 
eroclinic orbits or separatrices. One of them, which runs 
from {nd, 0) to (0, Trd) is shown in Fig. ^. It has the 
following time dependence. 
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the other "time surfaces" S* (see Fig. g) have a simi- 
lar functional form and an analogous dependence on the 
crossing times t„'s. 




where to is the moment of crossing the surface of section 
Eq (see Fig. 0). The heteroclinic orbits passing through 



FIG. 3. A sketch illustrating the construction of the 
separatrix map. Two unperturbed orbits, poit) and Pe(t) 
are shown. They follow two constant energy contours, 
[/o = (the separatrix) and Uq = e < 0, respectively. 
The energy-time coordinates e„ and f„ are defined by the 
crossings of the trajectories with the surfaces of section 
Eg and S^ (shown by bold segments). 

As explained in the Introduction, the unperturbed sep- 
aratrix is dressed with a narrow stochastic layer. In the 
case of the chessboard potential, this layer has a topology 
of a square network. We are interested in the long-time 
asymptotic behavior of the chaotic transport along this 
network. An efficifint tool to study such a transport is the 
separatrix mapE^'E^I The separatrix map is an approxi- 
mate map describing the dynamics near the separatrix. 
The application of the separatrix map to transport prob- 
lems has been previously considered in Refs. ES-E^. 

To construct the separatrix map we will consider "en- 
ergy surfaces" S^ in addition to the introduced above 
time surfaces E* . To avoid confusion we will elaborate a 
bit on the definition of such surfaces. E^'s and E* 's are 
introduced for each chessboard cell. Index q runs from 
to 3. The energy surfaces come through the saddle points 
and the time surfaces are drawn through the links con- 
necting the neighboring saddle points. The locations of 
E^'s and E*'s near the perimeter of the cell at the ori- 
gin are clear from Fig. S. The locations of the surfaces 
of section in the other cells can be obtained by periodic 
translations. Thus, index q in E* refers to the position 
of the corresponding link with respect to a given cell's 
center. Similarly, index q in E^ refers to the position of 
the saddle point. 

Let p{t) be the exact trajectory near the separatrix. 



As t increases, p{t) crosses the surfaces S^ in certain or- 
der. We denote by g„ the index of S^ at ri-th crossing and 
by e„ the value of Uq at this moment. Due to the time- 
dependent terms in the Hamiltonian, e„ changes with n. 
Let us find the difference e„+i — e„. The time derivative 
of Uo is given by 



dUo _ 1 
dt rriLO, 



E 

fc#0 



dUg dUk dUo dUk 
dpy dpx dpx dpy 



-ikuJct 



Ah f/fc's in this equation have to be calculated on the 
exact trajector y p{t) , which is not known. Therefore, 
following Refs. 0,pl^,p2l we perform the following ap- 
proximations. First we replace the exact trajectory by 
the unperturbed one with Uq — e„. Second, having in 
mind that |e„| <C Wq, we replace the trajectory with 
Ua = Cri by the separatrix motion po (t — t n) where pa{i) 



is given by equations similar to Eq. (2.11) and i„ is the 
moment of time when pit) crosses the surface of section 



S* . As a result, we find 



Cn+l 



where M„ is given by 



M„(i„), 



(2.12) 



M„(i) 



-E 



mojc 



k^o_ 



,.n dUodUk_ dUodUk\ ,ku^t' 

dpy Opx OPx dpy 



and is termed the Melnikov functionHS The integration 
is done along the trajectory po(i' — t). It can be shown 
that the terms fc = ±1 in the sum yield the dominant 
contribution. Thus, the Melnikov function can be ap- 
proximated by the integral 



M.it)^2,.MRefdt' '^^^;^-^]'-'^, 
J cosh[7u;c(i' - t)] 

— OQ 

X JiiRc/d) exp(^-iiu^t' + J + ^ 



(2.13) 



The integral can be evaluated by shifting the integra- 
tion path to the complex plane of t. Then Mn{t) can 
be represented by the sum of residues at the poles of the 
integrand. The residues from the poles closest to the real 
axis dominate the sum. Retaining only these terms, we 
arrive at 



M„{t) ~ Aesiwdn, 



dr, 



MRc/d) 
IT Trq„ 



i^ctn + -r + 



(2.14) 
(2.15) 

(2.16) 



Combining formulas (2.12) and ( ^.14 ) we obtain the first 
equation of the separatrix mapping 



e„+i = e„ + Aesin?9„. 



(2.17) 



To have the mapping in a closed form we need another 
equation relating i„+i to t„ and e„. Following Refs. 



y,ml,Ea, we take 



tn+l ~ tn + T-^(^n+l)' 



(2.18) 



where T(e) is the period of the unperturbed orbit 
Uq{p) — e. K straightforward computation gives 



T{e) 



1 

1 



-K 1 






In 



8Wn 



\e\ « Wo, 



(2.19) 



K being the complete elliptic integral oLt 

ceBE 



irst kind. 
Although it is a common practicdlitEiTEj to make 
the approximations similar to those we made above, 
their validity is far from being obvious. The justifica- 
tion has capie only recently with a new development by 
Tresche-tEa The extension of Treschev's analysis to our 
problemEJ indicates that the naive calculation of the Mel- 
nikov function is correct for Re <C d. If Re ^ d, then 



Eq. ( p.l5[ ) is off by a numerical factor and the replace- 
ment 



JijRc/d) 
JoiRc/d) 



jiRc/d) 



(2.20) 



is needed, where j{x) is some function of the order of 
unity for all real x. 

Besides the analytical methods, the validity of the sep- 
aratrix pian has been investigated numerically by several 
authorsESEIl and has been rated from "satisfactory" to 
"excellent." In the rest of this subsection we will assume 
that this is the case and calculate two quantities relevant 
for the transport, the width Acwcb of the stochastic layer 
around the separatrix and the average diffusion coeffi- 
cient D. 

We estimate Ae„ob following Ref. 02. First, we note 
that the relative change in e„ after one application of 
the separatrix map is small provided |e„| ^ Ae. Under 
this condition the map can be linearized. The defining 
parameter of the linearized map is /C, 



/C = 



1 



fd^n, 



COa-dn V ^'^n 



-1 = 



o^cAe dT{en) 
4 dcn 



Ae 

7en' 



and the map itself coincides with the standard map.EHI 
The crossover to the global stochasticity in the standard 
map occurs at |/C| ~ 0.97 [Ref. |31[], which yields the esti- 
mate 



eweb ^ I8j [Re/ d) 



2j2 



mujed 

7 



-ir/27 



(2.21) 



for the stochastic layer's width. Note that Cwcb ^ ^^/l 
is much larger than Ae, and so the approximation by the 
standard map is justified. 



Let us now turn to the evaluation of the diffusion coef- 
ficient D. For the chessboard geometry this wobleni has 
been considered previously by Ahn and Kim.EZi Unfortu- 
nately, they calculated the diffusion coefficient averaged 
only over the trajectories inside the stochastic layer. We, 
however, are interested in the diffusion coefhcient aver- 
aged over the entire phase space. Our approach to cal- 
culating D is close in spirit to the ones used for calcula- 
tion of the diffusion coefficient in planar periodic vortical 
flows, e.g., Rayleigh-Benard cells.LJO The details of the 
calculation can be found in Appendix H. The result is 



D = 0.45- 



mu>c 



Substituting this value into Eq. (2.22), we obtain 



D = 7.9 j {Rc/d) Lo^cf e-'^l'^'^. 



(2.22) 



(2.23) 



Note that apart from the numerical factor, D can be 
obtained from following simple arguments. Consider an 
ensemble of particles moving in the chessboard poten- 
tial. Their diffusive motion can be visualized as a random 
walk from one chessboard cell to the next. The motion of 
each particle is a combination of the drift along the cell 
perimeter and the series of random displacements in the 
transverse direction. The rate of diffusion depends on 
the distance of a particle from the cell boundaries. The 
particles located within a distance of one transverse step 
from the cell boundaries possess the fastest rate because 
they can cross to the neighboring cell after a single pas- 
sage along the cell's side. Particles further away from the 
perimeter remain trapped within the same cell for much 
longer time. Hence, their diffusion rate is negligible. Nat- 
urally, we can consider a model with an e-dependent dif- 
fusion coefficient D{e) = e(Ae- |e|) d§/r(e) where Q){x) 
is the step- function and do = V2T:d is the length of the 
cell's side. The net diffusion coefficient can be obtained 
by averaging D(e) over the phase space, i.e., over the area 
in coordinates {px,Py), 



D^ydeDie) 



dSje) 



where S{e) is the area of the cell's region bounded by the 
contours Uq = and Uq = e. It is trivial to show that 
dS{e)/de — T {e)Jjnuj c] therefore, D — Ae/mtjj,: which 



reproduces Eq. (|2.22| ) up to a numerical factor. 

Finally, the diffusion coefficient can be written as a 
function of the magnetic field B, 



In 



D 



to. 



B 

'b~^ 



3/2 



(2.24) 



where 



2^6 V^^I^^ /M^y/^ 



[cf. Eq. (pT7|)]. Formula (2.24) was derived assuming that 
7 <C 1, i.e., that B ^ B^h- In addition, we assumed that 
Re > d, which is equivalent to S < B^h {E/W)^^'^. As 
one can see, the dependence of 13 on i? for the chessboard 
geometry is given by a squeezed exponential with the ex- 
ponent I . In the next subsection we treat a more general 
case of a Gaussian random potential. We will show that 
the squeezed expo nential is replaced by a simple one as 
given by Eq. (f 



B. Gaussian random potential 

A Gaussian random potential is fully specified by its 
two-point correlator C(ri — r2), 

C(ri-r2) = (C/(ri)C/(r2)), 

C(0) EE W^. 

In many cases, it is also convenient to deal with the 
Fourier transforms of C/, which have the following cor- 
relator 

{U{qi)U{q2)) = {27r)H{q,+q2)C{qi) 

(Fourier transforms are denoted by tildes). Given the 
function C{r), we want to calculate the diffusion coef- 
ficient in strong magnetic fields. Similar to the case of 
the chessboard potential, let us first investigate the "un- 
perturbed" motion, the drift along the contours Uq{p) — 
const. Clearly, Uo{p) is also a Gaussian random potential 
with correlator Co related to C by 

Co{q)^[MqRc)fC{q). 

The unperturbed motion is determined by the prop- 
erties of the level lines of Uq. It is known that all of 
such lines except one, the percolating contour, are closed 
loops. The Gaussian random potential shares this prop- 
erty with the chessboard potential considered above. In 
addition, the position of the percolation level is the same 
for both potentials: Uq = 0. There exists, however, an 
important difference in the properties of level lines in 
the two cases. The diameters of the loops in the chess- 
board do not exceed 2TTd. On the other hand, constant 
energy contours of the random potential can have arbi- 
trarily large diameters. Such large loops are found in the 
vicinity of the percolating contour. (The latter one can 
be considered as a loop with infinitely large diameter). 
As the diameter of the contour increases, the range of Uq 
found at such contours shrinks, tending to the percola- 
tion level Uq = 0. 

Similar to the chessboard geometry case, the exact tra- 
jectories do not simply follow the level lines of Uq{p) but 
exhibit small transverse deviations from them. As a re- 
sult, a finite diffusion coefficient appears. To calculate D 
we will use a close analogy of the problem at hand with 



the problem of calculating the effective diffusion coru 
stant of a particle diffusing in an incompressible flow.Efl 
Below we essentially reproduce the basic arguments of 
Isichenko et alE3 with slight modifications appropriate 
for our problem. 

Borrowing the terminology of Ref. |3J, we call a bundle 
of constant Uq contours with diameters between a and 2a 
a convection cell or an a-cell. The values of Uq in typical 
a-cells belong to an interval [— w(a), w{a)], which narrows 
with increasing a. Let us denote by L{a) the perimeter 
length of typical o-cells and by Ae(a) the change in Uq 
accumulated along the trajectory following the perime- 
ter, for which the time T{a) '^ L{a)/vd is required. The 
key point in estimating I? is a ramification between mix- 
ing [with Ae(a) > w{a)] and non-mixing [Ae(a) < w{a)] 
cells. It takes a single period T{a) or even a fraction of 
thereof for the particle to leave a mixing cell, whereas 
particles in non-mixing cells remain trapped for time in- 
tervals much larger than T{a). The dominant contribu- 
tion to the transport comes from the mixing cells of the 
largest width w{a) for which Ae(a) ~ w{a). Denote the 
diameter by such cells by Um- The particles situated in 
such cells perform a random walk from one optimal cell 
to the next. The characteristic step of the random walk 
is Qm and the characteristic rate of the steps is 1/T{am)- 
Thus, the diffusion coefficient of such "active" particles 
is of the order of af^/T{am)- The net diffusion coeffi- 
cient can be found by multiplying this diffusion coeffi- 
cient by the fraction of the total area occupied by the 
optimal convection cells. Note that the width of the a-- 
cells in the real space is of the order of Ae(am)d/Wo. 
Using this, the fraction of the area can be estimated to 
be [Ae{am)d/WQ]L{am)/al^ = [Ae{am)/muJc]Tiam)/al,. 
Finally, we obtain 



D 



Ae{ar,i) 



mojc 



(2.25) 



which closely resembles Eq. ( 2.22 ) for the chessboard.til 
However, now Acm, = Ae{am) depends on the diameter 
a,n of the optimal cells, which has yet to be found. We 
see that the calculation of D hinges upon the calculation 
of Aem • To accomplish the latter task we can make the 
same kind of approximations as in deriving the separa- 
trix mapping for the chessboard. Then we obtain the 
following expression. 



Ad = E 1^ 



2 



riT^O 



A„ EE (bdtVd'VUn[po{t)]e' 



(2.26) 



(2.27) 



where the integration path is the unperturbed orbit 
Uo[po{t)] = const belonging to a given Om-cell. Observe 
that the integrand is the product of a slowly changing 
function /„(i) = Vd^Un[po{i)\ and a rapidly oscillating 
exponential factor e^™'^'=*. It is customary to estimate 
such integrals by shifting the integration path into the 



lower half-plane of complex t where the oscillating factor 
decays exponentially. By using the method, one arrives 
at the following estimate 



Ae 



|Ai| 



fe 



^27rii?fee-|i™^''l'^'^ 



(2.28) 



where r^ are the singular points of the function /i(t) 
in the lower half-plane plane and Rk are some pre- 
exponential factors. For example, if /i(i) has a simple 
pole at Tfc , then R^ is up to a phase factor the residue of 
such a pole. Equation ( 2.28| ) is similar to Eqs. (2.14-2.16) 
for the chessboard potential. 

Denote the coordinate along the drift trajectory by s, 
then /i(t) = Vd (dUi/ds). The singularities of fi{t) may 
originate either from Vd or from (dUi/ds). Let us inves- 
tigate the former possibility. To get the necessary insight 
we will use the exactly solvable model of the chessboard 
potential, which we studied above. In the latter case 



Vd{t) - 



\/2juJcd 



cosh[jujcit~ to)] 



(2.29) 



[see Eq. ( 2.1l| )] and the singularities of Vd{t) in the lower 
half-plane consist of the "parent" pole at tQ — iTT/2^LUi. and 
a series of "daughter" poles at ip — itt [k + i) /■fWc, k = 
1,2.... Note that the imaginary part of the parent pole 
is of the order of the characteristic time scale {-^ojc)^^ of 
the drift motion. 

In the case of the random potential, we also expect to 
find series of singularities of Vd{t). However, there will 
be not a single series but a large number N{am) of them. 
Indeed, Vd{t) has about L(am)/d minima on the trajec- 
tory s{t) . The points of minima divide the trajectory into 
L{am)/d intervals of length ~ d. In each interval Vd{t) 
first rises, then reaches a maximum, then decreases, i.e., 
it exhibits the same kind of behavior as in the chess- 
board case. Therefore, a naive estimate of Niom) is 
N{am) ^ L{am)/d. Since ImTk^s enter Eq. ( 2.28| ) in 
the arguments of the exponentials, the dominant contri- 
bution to Aem comes from these N(am) parent singu- 
larities. Let us now discuss Imr^'s. It is obvious that 
different am-cells give rise to different Imr^'s, i.e., there 
exists a certain distribution of Imrfc's. What kind of dis- 
tribution should we expect? Clearly, the typical value of 
the imaginary parts of the parent singular points should 
be of the order of the characteristic time scale of the drift 
motion, {'yujc)~^, where 7 can be defined as follows: 



Wn 



1 



mui, 



H^' 



with Wo and d being 



Wo = V^o(O), d 



Co 



2V^Co' 



However, it would be a mistake to think that Aem is de- 
termined by this typical value. Indeed, the deviations 



of Imrfc from their average value are dramatically en- 
hanced in Aem- Therefore, we can expect an extremely 
broad range of the exponentia l fact ors entering the sum 
on the right-hand side of Eq. ( |2.2^ ). At the same time, 
there is no such an enhancement for R^. This kind of 
arguments imply that we can estimate Ae^ considering 
only the distribution of Imrfc 's, i.e.. 



Aei 



^e*>?.g-|i. 



Tfcll^ 



where dk is the phase of the complex number R^. Of 
course, the pre-exponential factor in Acm can not be 
found by this approach. We will further assume that 
i^fc's are uncorrelated, which results in 



Aei 



^e-li. 



Tfcl" 



From this, we find that 



Ae 



W^ 



2L(a„ 



d7'P(7')e-2/^', 



(2.30) 



where 7' = l/ImrWc and P{^') is the distribution func- 
tion of 7'. The fist factor on the right-hand side is written 
solely to provide the correct dimensionality. 

It is possible to show that P(7') has the Gaussian tail. 



P(.l') 



exp 



AY' 



7'>7, 



(2.31) 



where A '--^ 1 is some number. This result can be obtained 
from the following simple physical arguments. More rig- 
orous treatment is relegated to Appendix g. 

Let us again look back at the chessboard model. As 
one can see from Eq. ( 2.29 ), Vd as a function oft exhibits 
a brief pronounced pulse near its maximum at t = to- 
The duration of the pulse is of the order of {jujc)"^- 
It is this time scale that determines the imaginary part 
of the closest singular point. Let us now return to the 
random potential case. One can speculate that singular 
points of Vd{t) are always associated with such kind of 
pulses. By this argument, the singularity at the point 
ts = ti — it2 with < ^2 <^ il^c)~^ requires an un- 
usually short pulse of duration At ^ ^2- To produce 
such a pulse Vd{s) must have a large and sharp maxi- 
mum. Let us estimate, e.g., the height of this maximum. 
The half-width As of the maximum is of the order of 
As ~ \j^^dMd- O^ the other hand, we should have 
As ~ Vdt2- Thus, Vdv'lj^ ~ —t^ , which shows that small 
values of Im t^ require large values of Vd and its second 
derivative, Vd ^ djti and w^' ~ l/t2d. Recall now that 
the distribution functions of both Vd and v'^ have Gaus- 
sian tails, so that the probability of finding an unusually 
large Vd is of the order of exp(— Aiw^/7^(jj^(i^) and sim- 
ilarly for v'^ [Ai ~ 1 is some number). Substituting 



d/t2 ^ j'iOcd for Vd, we arrive at Eq. ( 2.31 ). The calcu- 
lation of A for some particular example of C{r) can be 
found in Appendix H. 

The estimation of the integral in Eq. ( 2.30] ) by the 
saddle-point method results in 



Aei 



Wr 



L(a„ 



■ exp 



3^1^ 

" -72/3 



(2.32) 



On the other hand, L(am) obeys the scaling law 

L(a„)(x|Ae„r'"^\ (2.33) 

where v and dh are some exponents, which depend on 
the properties of the correlator Co{q) [Ref. Q. Their ac- 
tual values are not very important at this point. Equa- 
tions ( p. 32 ) and (2.33) enable one t o find Ae^, which can 
then be substituted into Eq. (2.25). As a result, we find 
the diffusion coefficient. 



D -- cucd'^'-f"' exp 



J 



7 



2/3 



(2.34) 



where a is some number and 

3^1/3 



J = 



1 + i^dfi 

is another number. We remind the reader that we can- 
not calculate the correct pre-exponential factor in for- 
mula ( 2.34 ). The particular choice of this factor made 
in Eq. ( ^.34| ) provides a matching of this equation with 
Drude-Lorentz formula (^^) at 7 = 1 where the both 
formulas give D ~ cocd'^ (up to p ur ely numeri cal fa ctors) . 

|) if one 
for the 



1.6), and 



This can be seen from Eqs. ( |l.2| ), 

takes into account the approximate expressior 

transport time t. 



E 



d 



In this subsection we implicitly assumed that t he i n 

3/2 



equality Rc'3> d holds. In t his ca se ^ oc B '^/•^ [Eq, 
Substituting this into Eq. (2.34), we obtain 



m 



D(xucd'^-f"e'^/^'', B > Be. 

declared previously in Sec. ^ Note that the dependence 
of D on the magnetic field is given by a simple expo- 
nenti al no t the squeezed one as in the chessboard model 
[Eq. (2.24)]. The reason for this difference comes from 
the important role of rare places on the trajectories with 
unusually sharp features of the averaged potential Uq. 
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III. INTER-LANDAU LEVEL TRANSITION 
AMPLITUDES 

In the preceeding Section we showed that in strong 
magnetic fields, B > Be, the guiding center of the cy- 
clotron orbit closely follows the level lines Uq — const 
of the averaged potential Uq. Nonvanishing diffusion co- 
efficient appears due to small deviations from the level 
lines. The characteristic value Acm of such a deviation 
was calculated purely classically. Due to the energy con- 
servation, Ae„i also represents the change in the kinetic 
energy lujc of the particle [Eq. (2.2)]. 

The purpose of this Section is to calculate the change 
in kinetic energy quantum-mechanically by taking into 
account the discreetness of the spectrum, i.e., the exis- 
tence of the Landau levels (LLs). Note that this is not 
yet a consistent quantum-mechanical treatment of the 
problem. For example, in this Section we ignore localiza- 
tion and/or quantum tunneling. An attempt to touch on 
some of those complicated issues will be postponed till 
the next Section. 

In quantum- mechanical terms, the change in kinetic 
energy results from intcr-LL transitions. Indeed, the 
change in kinetic energy due to N —> N + k transition 
is equal to ktiiOc- Denote the transition amplitude upon 
the completion of the loop Uq ~ const by Apf^pj^^, then 
(^Cm) is given by 



(Ad 



{huJc)^'^k'^\AN,N+k\ 



(3.1) 



It is obvious from this formula that the inter-LL tran- 
sitions may be significant only within a certain band of 
LLs. If Acm is larger than huc, then the number of LLs 
in that band should be of the order of Acm/^ojc- Denote 
by _B* the field where Ae™ = fiujc- In fact , this notation 
has already been used in Sec. | [Eq. ( |1.10| )]. If S > S*, 
then Ae,„ < hujc and even the transitions to the neigh- 
boring LLs must be suppressed. In this case the sum over 
k is dominated by the two terms, k = ±1; therefore, 



1^ 



(Ael 



N,N±1\ 



2{nuocf 



(3.2) 



In deriving Eqs. (3.1) and ( |3.2[ ) we implicitly assumed 
that the classical and the quantum calculations of (Ae^) 
give the same result. This will be demonstrated below. 
Before w e do so, let us mention one interesting fact. 



Using Eq. ( 2.25| ) and the Einstein relation axx = e^vD 



where u = mj'n'h is the density of states (de Haas- 
van Alphen oscillations neglected) one arrives at the fol- 
lowing formula 



Ae 



h hu)c 



It can be interpreted in the following way: the trans- 
port is determined by the aforementioned band of about 



(Aem/?iWc) LLs with energies near the Fermi energy. 

Each level contributes ^ to axx [cf. Refs. pTf . 

The general formula for A^^pf^f^ derived in Ap- 
pendix O reads 



A 



N,N+k — 



2^ 



e-*'=^exp 



Y^ An ~inB 



n=^0 



ihiJr 



(3.3) 



where A„'s are given b y E q. ( 2.27 ). Substituting this ex- 
pression into formula (3.1) and taking advantage of the 
identity 



2ir 



/^ E fcV^V(^) = -/"(o), 



k^ — oo 



we recover the classical formula ( ^.26 [ ) fo r Ae^. 

Finally, it is easy to see that Eq. (|3.2| ), which we de- 
rived without any calculations, is consistent with for- 
mula (3.3). Indeed, |Ai| ~ Ae^. If Ae„i <S^ hcuc, then 



the second exponential in Eq. (3.3) can be ex pan ded in 
the Taylor series, which trivially leads to Eq. (3.2). 



IV. 



CALCULATION OF THE QUANTUM 
LOCALIZATION LENGTH 



In Sec. we argued that the localization length is ex- 
ponentially large in weak magnetic fields and has to de- 
cay as the magnetic field increases. This statement is an 
oversimplification in two respects. Firstly, ^ is, in fact, 
expected to diverge at certain discreet values B^ of the 
magnetic field 



? = ^o 



B 



N+l 



B 



N 



B- B 



N 



(4.1) 



where fj, is a. critical exponent .u Secondly, such diver- 
gences neglected, £, starts decreasing only from B ^ 
hc/eltr, at which the magnetic length I = ^Jh/rruJc be- 
comes of of the order of the transport length /tr = vt. 

Let us discuss these issues in some detail. Scaling the- 
ory of localizatiop-is one possible way to approach this 
difficult problem.lia In scaling theory one tries to under- 
stand the localization by considering the behavior of the 
dimensionless conductance g = -^axx as a function of 
system size L. This behavior is described by the scaling 
function 



Pig) 



d\n.g 



(4.2) 



One starts with calculating the conductance go = 5(^0) at 
some short length scale L = Iq where it is large and then 
finds how g is renormalized towards larger L. The lo- 
calization length is the length scale where g{L) becomes 
of the order of unity. (If go is of the order of unity of 
smaller, then a different approach has to be used, see 
below) . 
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It has been conjecturedES that all physical system can 
be grouped into certain universality classes with the same 
functional form of the scaling function. If we neglect 
the spin-orbit coupling, then the appropriate universal- 
ity class for our system is determined by the relation 
between L and the magnetic length I. For L ^ I, the 
system belongs to the adipgonal class, where the scal- 
ing function is given byEj'Lj 



Co =Zoexp(7r 5o), go > 1 



(4.8) 



(3{g)^ , L«/. 

ng 



(4.3) 



For L ^ I the system is in the unitary class. The scaling 
function is given by 



Pia) 



1 



2n^9 



2„2 ■ 



L>L 



(4.4) 



The latter result was derived both by the conventiortal 
diagram techniqucc3c2l and by an effective field theory.E3 
Solving the scaling equation (4.2) for g{L), we find that 
C experiences a growth from the value of 



C 



Ztrexp(^- 



kpl 



FHi 



at B = to 



^-k 



. exp (w'^kpl^^) 



(4.5) 



(4.6) 



at _B ~ hc/el^j. where I = hr- In stronger fields, 
B > fic/ell^, the system belongs to the unitary class fXf^ 
relevant length scales and ^ is given by the formulallSliS 



C = Zoexp[^2go(5)'] 



(4.7) 



following from Eq. (4.4). The dimensionless conductance 
go{B) decreases with B. For the case of a long-range 
random potential this follows from the results of the pre- 
ceeding Sections. Therefore, the initial growth of f at 
very weak magnetic fields is followed by the exponential 
decay of ^ as -B increases. This is the statement we put 
forward in Sec. |. 

Unfortunately, Eq. (4.7) cannot be entirely correct 



beca use it does notj-ijeproduce the critical divergences 
[Eq. ( [4.l|) ]. PruiskerO argued that the critical behavior 
is a non-perturbative effect. His field-theoretical treat- 
ment yields an expression for the /3-function, in principle, 
different from the simple form (4.4). However, the devi- 



ations from Eq. (4.4) become significant only when the 
renormalized value of g app roaches unity. On this basis 
we speculate that Eq. ( ff.T] ) gives only the lower bound 
for the localization length ^. We further assume that this 
lower bound is close to the actual value of ^ away from 
criticality. In other words, Eq. ( [4.7| ) gives, in fact, not ^ 
itself but its non-critical prefactor ^o entering Eq. (4.1). 
Note that ^ ~ ^o at the midpoints between neighboring 
divergences of C, i.e., at the QHE conductivity minima. 
This is exactly the quantity discussed in Sec. m where we 
postulated the ansatz (1.14). We will rewrite it here for 
the ease of reading: 



By virtue of this ansatz^ the calculation of ^o boils down 
to the evaluation of the short length scale conductance 
50- 



Previous attemptstjcil to treat the localization prob- 
lem in the QHE have been focused on the case of a short- 
range random potential, i.e., the potential whose corre- 
lation length is much smaller than de Broglie wavelength 
27: /kp. In this case go has to be calculated quantum- 
mechanically, e.g.., wit hin a Self-Consistent Born Approx- 
imation (SCBA) .113113 Recall that our the ory a pplies to 
the case kpd > (E/W)^^^ > 1, see Eq. (|l.ll| ). There- 
fore, there is a whole intermediate region 1 <^ kpd <^ 
(E/W)'^^^ separating the domains of applicability of our 
and the previous theories. The calculation of ^o in that 
region is a separate problem and will be discussed else- 
where. 

In the case of long-range random potential, which we 
consider here, go{B) can be calculated with the help of 
Einstein relation, 



goiB) = hv{B)D{B), 



where v{B) is the density of states at the Fermi level 
and D{B) is the classical diffusion coefficient. Accord- 
ing to the results of the previous Sections, D{B) is given 
by Drude-Lorentz formula (1.2) at B < Be and by for- 
mula ( 2.34| ) at B > Be- Let us now discuss the behavior 
of I'iB). In principle, i^iB) oscillates with B around its 



zero field value iy{0) — m/irfi . However, for B smaller 
or at least not to much larger than B^ such oscillations 
are exponentially small because the width of LLs, which 
is of the order of Wo [Ref . M , is much larger than the 
distance fiUc between them. Therefore, we can use the 
zero field value v{Q). 



Substituting all these results into Eq. (4.8), we obtain 
£,o{B). The f unctio nal f orm of this dependence is given 
by Eqs. ( |1.15 ) and (1.16). Graphically, it is illustrated by 
Fig. [l[ Observe that the overall decay of ^o as B increases 
becomes extremely sharp at B > Be- As a consequence, 
already at the field B — B^, , which is only logarithmically 



larger than Be [see Eq. ( |l.7|) ] ^o ceases to be exponentially 
large. At B > B^ go becomes less than one and Eq. ( |4.8| ) 
does not hold any more. In this region the localization 
length is determined mainly by quantum tunneling rather 
than by the destructive interference of classical diffusion 
paths. Thus, the calculation of ^o requires a different 
approach. It will be discussed in a forthcoming paper to- 
gether with the prefactor in formula (1.16). At this point 
we can only say that ^o is expected to have a power-law 
dependence on B and eventually match the predictions 
of Raikh and ShahbazyanEl at sufficiently large B. 
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V. DISCUSSION AND CONCLUSIONS 

In this paper we studied a two-dimensional motion of 
a charged particle in a weak long-range random poten- 
tial and a perpendicular magnetic field. We showed that 
the phase-space averaged diffusion coefficient is given by 
Drude-Lorentz formula only at magnetic fields B smaller 
than certain value Be- At larger fields, the chaotic mo- 
tion is suppressed and the diffusion coefficient becomes 
exponentially small. 

To make connection with the experiment our results 
can be applied to the following model. We suppose that 
the random potential is created by randomly positioned 
ionized donors with two-dimensional density rii set back 
from the two-dimensional electron gas by an undoped 
layer of width d. We will assume that mtP 3> 1 and 
also that d^ as, where a^ is the effective Bohr radius. 
In this case the random potential can be considered a 
Gaussian random potential whose correlator is given in 
Appendix |Q. As a particular example, we consider a spe- 
cial case where the density of randomly positioned donors 
is equal to the density kp/{2TT) of the electrons. We call 
it the standard potential. It is easy to see that for the 
standard potential E/W ~ kpd and the domain of appli- 
cability of our theory [Eq. (1. 11)] is simply kpd ^ I. In 
modern high-mobility GaAs devices this parameter can 
be as large as ten. It is easy to verify that the mag- 
netic field Be where the classical localization takes place 
corresponds to LL index Nc ~ {kpd)^/^^ which can be 
a number between 10 and say, 50 for GaAs heterostruc- 
tures. Another important magnetic field B* [Eq. (1.10)] 
corresponds to LL index A^* , which is only slightly smaller 
than Nc- As explained in Sec. ||, Af, is the number of the 
"first" QHE plateau in the sense that observability of 
plateaus with larger A^ require exponentially small tem- 
peratures. 

The point A^ ~ N^ plays another important role. It 
is the largest A^ where it is possible to see the activated 
transport o'xx oc e~^"'"^, Ea — fiujc/'^ in the minima 
of a XX- Indeed, it is known that in strong fields or for 
small A^'s the dissipative conductivity demonstrates the 
Arrenius-type behavior at not too low temperatures. As 
the temperature decreases, the activation becomes re- 
placed by the variable-range hopping, see Eq. (1.12). 

Equating the two exponentials, we find the tempera- 
ture Th at which the activation gives in to the hopping. 



Th 



{flLOc 



To 



This formula can also be written in another form. 



huJr 



— const , 

To TsRc 



(5.1) 



(5.2) 



where Vg — \f2e^ JKhvp is the gas parameter, which is 
of the order of unity in practice. Let us demonstrate 
that the activated behavior should not be observable at 



B < B^- Indeed, it make sense to talk about the acti- 
vated behavior only at temperatures below the activation 
energy Ea — huoc/'i- Therefore, the activated transport 



can be observable only if the right-hand side of Eq. (5.2) 
is less than unity. Thus, in the vicinity of B — B^ or near 
A^*-th conductivity minimum where ^q suddenly becomes 
exponentially large with decrease in B, there is no place 
for the activated transport left. Strictly speaking, this 
argument only proves that the activated behavior is ab- 
sent at i3 ^ B*. However, it can be shown, and it is a 
subject of a forthcoming paper, that the ratio £,q{Bc)/ Re 
is smaller than one in the standard case. Therefore, the 
activation indeed disappears at -B ~ i?* rather than at 
some much stronger field. 

The behavior of ^o in magnetic fields stronger than B^, 
has not been investigated in the present paper. It will be 
discussed elsewhere. We expect that at such magnetic 
fields (.o{B) is a certain power law matching the results 
of Raikh and ShahbazyanEd at sufficiently large B- As 
explained in Sec. ||, such a dependence can be studied 
experimentally. 

Finally, in this paper we have neglected the influence 
of electron-electron interaction on ^o • This complicated 
issue warrants further study. 
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APPENDIX A: DIFFUSION COEFFICIENT IN 
THE CHESSBOARD MODEL 



To calculate the numerical factor in Eq. ( 2.22 ) for the 



diffusion coefficient we proceed as follows. First, we will 
introduce the random-phase modecj arguing as follows. 
The well-known property of the standard map is a fast 
mixing in the phase variable z?. The correlations in phase 
decay according to (e^(''"-''«)) ~ |A:|-"/2 [Ref. O as a 
function of the iteration number n; therefore, for \K\ ^ 1 
the phase memory is typically lost after a single itera- 
tion of the map. The situation with the separatrix map 
is similar, which allows a simplification of the problem. 
We will a ssume that e„ is still transformed according to 
Eq. ( 2.17 ) as long as |e„+i| < ewcb- If the new value of 
|e„+i| is larger than Cwcb, then e„+i = e„. At the same 
time, "dn will be a purely random variable uniformly dis- 
tributed in the interval (0, 27r). As we will see below, 
for transport only the narrow boundary layer |e| '^ Ae 
is important (cf. Refs. ^|3|) where jATj :^ 1 and there- 
fore, such a random-phase model is adequate. Ahn and 
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KimEZi studied this model numerically and found an ex- 
cellent agreement between the diffusion coefficients found 
from the random-phase model and from the original sep- 
aratrix map. (Of course, the random-phase model lacks 
certain features of the original separatrix mapping, e.g., 
a rich hierarchical island structure). 

Consider now an ensemble of particles, each having 
the same total energy E but different initial conditions 
at i = 0. In the original problem with Hamiltonian (2.7), 
we can describe this ensemble by a distribution function 
(guiding center density) F{p,t). We will calculate the 
diffusion coefficient as the coefficient of proportionality 
between the the average particle flux and the average 
gradient of F in the stationary state. It is convenient 
to rotate the coordinate system by j. We denote new 
coordinates hy £ and rj. The gradient of F is in the r)- 
direction (Fig. m. 



faces of section S^. The superscripts distinguish between 
the positive and negative e contours (Fig. Q). Func- 
tions F^ are taken to be zero for e < and similarly, 
F~ (e) = for e > 0. We also define "full" functions F„ 
byF„(e) = F+(6) + i^-(e). 

Denote the length of the chessboard cells by do {do = 
^/2-Kd) and the average gradient (|Vi^|) by AF/do, then 



D 



$ $1 - $0 

AF ^ AF '■ 



(Al) 



where $ is the total flux through the side {0 < ^ < 
do, 77 = 0} and $„ is the flux incident upon the S^ 
surface, 



$„ 



l Wo 

dl vail) F+[e{l)] = ^ f deF+{e). (A2) 
muJc J 
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FIG. 4. The chessboard in the coordinate system ro- 
tated by .|. Pluses and minuses at the centers of the 
chessboard cells label maxima and minima of the poten- 
tial. The direction of the drift velocity is indicated by 
arrows. Bold segments are the surfaces of section S^ , the 
same as in Fig. H. Distribution functions F^ represent 
the deviation of the guiding center density from the sam- 
ple averaged value at those parts of S^'s, which are inside 
of the two cells on the right. Surfaces Eq and Ef also 
penetrate the two cells on the left. The corresponding 
distribution functions are related to Fj , F^ as shown. 



In fact, the description of the ensemble by function 
F, which is a function of a vector argument, is reason- 
able only when we study the exact dynamics. After we 
have replaced the exact dynamics with that of the sep- 
aratrix map and now even of the random-phase model, 
this kind of description became too detailed. Instead, it 
is sufficient to introduce a set of the distribution func- 
tions F^ie) of a single argument. Each function in the 
set represents the deviation of F from its average value 
at the intersections of the contour Uq = £ with the sur- 



Here I is the coordinate along S^ and Vd{l) is the drift 
velocity. 

To obtain the equation for F„'s note that within the 
random-phase model, Fn+i is quite simply related to F„. 
For example, 

F2{e) = J ^Fi(e - Aesin,?) - ^ o Fi(e), (A3) 



Similarly (see Fig. ^, 

AFe(-e) + F3-(e)+F+(e) 

= So[AFe{-e) + F,-ie) + F+{e)], (A4) 

where &{x) is the step-function. Suppose that all F„'s 
are equal to zero at the center of the cell, then the chess- 
board symmetry dictates F3 = — Fi and F2 = — Fq and 
also that functions F„ 's a re even. These relations can be 
substituted into Eq. (A4). Then one can eliminate Fq 
and obtain an equation solely for Fi , 



{1 + IoSoIo S')Fi(e) = (/ o 5 - /)AFe(-e), (A5) 

where (/ o /)(e) = sgn(e)/(e). Equation ( |A5| ) is an in- 
tegral equation, in principle solvable by the Winer-Hopf 
method. However, we have not been able to find its solu- 
tion analytically. At the same time, a numerical solution 
can be obtained rather easily. The result is shown in 
Fig. 1]. By calculating the area bounded from above by 
the graph of function Fi from below by the graph of Fq 
and from the left by the vertical line e = we ha ve ob - 
taincd the numerical factor 0.45 in the expression ( 2.22| ) 
for the diffusion coefhcient. 

Both Fo(e) and Fi(e) decay exponentially at |e| » Ae. 
This is in accordance with the statement above that 
only a narrow boundary layer is important for the trans- 
port. Sipulaj to the conventional advection-diffusion 
problems J^Ej the width of this layer, dAe/Wo, is of the 
order of the average displacement of the particle in the 
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direction perpendicular to the flow upon travelling the 
length of the chessboard's cell. 





^-0.2 



FIG. 5. The distribution functions Fq and Fi. The 
density (vertical axis) is in units of AF and the energy 
(horizontal axis) in units of Ae. 



APPENDIX B: ASYMPTOTIC BEHAVIOR OF 
FUNCTION P(7') 

In Sec. H we showed that the calculation of the diffu- 
sion coefficient in strong magnetic fields can be reduced 
to the calculation of the distribution f uncti on P{"i')- Our 
deri vation of the subsequent formula ( 2.34 ) was based on 
Eq. (2.31), which we rewrite here for convenience, 



Pii) ex exp ( ^ 



(Bl) 



Our goal is to derive Eq. (Bl) 



The sufficient conditions for Eq. (Bl) to hold seem to 
be as follows. Suppose that C{r) is isotropic, i.e., de- 
pends only on r = ^J x^ -|- j/^. Function C{r) must be 
analytic for all real r. In addition, C{r) must be analytic 
in some complex neighborhood of r = 0. Note that such 
conditions can be met only if C{q) decays s uffi cientl y fas t 
at large q. The "realistic" potential [Eqs. (CI) and (C2)] 
meets the aforementioned requirements. It is easy to see, 
for instance, that in this particular example C{r) is an- 
alytic within the circle |r| < 2(i in the complex plane of 
r. 

Below we use the symbol v[s), which is a concise nota- 
tion for the drift velocity on the appropriate closed-loop 
trajectory. The argument s has the meaning of the co- 
ordinate along this trajectory. Let us introduce the fol- 
lowing definition. Consider an arbitrary point r on the 
trajectory. We call the function u(s). 



u{s) 






n=0 



the maj orant of v{s) at the point r. The derivation of 
Eq. (Bl) is based on the following 



Lemma Suppose that function v{s) is analytic in a com- 
plex neighborhood of point r. Consider the solution s(t) 
of the Cauchy problem 



ds 



w(s), s{to) = r 



(B2) 



as a function of the complex argument t and define 
v(t) = v[s{t)]. In the similar way we define function 
u{t) where u{s) is the majorant of v{s) at point r. 

Under such conditions the shortest distance ty from ip 
to a singular point of v{t) is greater or equal to the sim- 
ilarly defined distance t„ for the function u{t). 



Corollary 1 i„ satisfies the inequality 



^M 



tv > 



ds 



Oq + ai{s - r) + 02(5 - r)2 + 



(B3) 



where 



L(»)(r)| 



and rm {r) is the smallest real number larger than r such 
that u{rm) = 00. 



Proof. The distance t^ can be calculated as 

l/n 



ty — lim inf 



n\vo 



ar-it) 



(B4) 



where the derivatives are taken at t = to and vq is an 
arbitrary positive parameter with dimension of velocity. 
A similar expression can be written for i„. Clearly, 



d^ 
dF- 



vit)U^ 



t=to 



ds 



yis)\s^r 



Upon taking the derivatives, the right-hand side becomes 
a sum of products of v^'''>{s). It is important that all the 
coefficients in the sum are positive, which leads to the 
inequality 





< 


d"u 


dt" 



The statement of the Lemma follows from this inequality 
and Eq. (g). 

To prove Corollary |l| note that the differential equation 
on s(t) with u{s) on the right-hand side can be solved in 
quadratures. 
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t{s) - to 



ds 

u{s) 



At singular points u(t) diverges; therefore, the distance 
to such a point can be found by choosing the upper hmit 
from the condition \u{s)\ = oo. It is easy to see that 
by choosing r^ for the upper hmit and the integration 
path along the real axis, we fi nd th e shortest distance i„ 
to the singularity. Inequality (B3) trivially follows from 
this result ■ 

Recall that we are interested in finding the imaginary 
parts of the singular points. For this purpose the follow- 
ing Corollary is helpful. 

Corollary 2 Suppose that v{s) is a periodic function 
with period L. The imaginary part Imtsof any singu- 
lar point tg of function v(t) defined by (B2) satisfies the 
inequality 



^ir)- 



Imts > min 

0<r<L 



ds 



flo + ais + a-is^ + 



(B5) 



The proof is quite obvious. 

At this point we are ready to derive Eq. (Bl). The 
derivation is based on the rigorous bound (B5) supple- 
mented by a few quite plausible assumptions. Consider 
some contour L'o(s) = const with perimeter length L and 
choo se a point s = r on this contour. All a„'s enter- 
ing ( [B5[ ) are random numbers with certain distributions 
such that a„ larger than some characteristic values A^ are 
exponentially rare. Generally, An is of the order of the 
typical value of (n -I- l)-th order derivative of Uo{p)/'muJc 
with respect to px and/or py, e.g., 9"+^C/o(p)/5p"^^. 
The latter quantity has the normal distribution with vari- 



ance 0{n ^) 



C, 



(2n+2) 



(0) 



Hence, 



C, 



(2n+2) 




(0) 



where a is some number. 

Let pm be the shortest distance from the origin to the 
singular point of Co(p) in the complex plane of p. W e 
will assume that pm < oo as in our example [Eq. (|C2| )] 
where pm = Sd. (The analysis of the other case, pm = cx3 
is similar). Based on this example, we further assume 
that Pm ~ d. 

Since pm is the radius of the analyticity circle of Co [p) , 
we must have [cf. Eq. (B4)] 



lim sup 



1 c(")(o) 



Co(0) 



l/r 



1 
Pm 



Consequently, the asymptotic behavior of An at large n 
is given by 



Ar. 



^LUcd 



1/n 



1 

d' 



(B6) 



The typical value of a„ in formula (B5) is of the order of 
An', hence, the typical distance of singular points of v{t) 
from the real axis is of the order of (tWc)"^- This is the 
statement we put forward in Sec. [IB. Next we would 
like to identify the condition for v{t) to have a singular 
point ts with |Imis| ^ (t^c)^^. It is more convenient to 
work with dimensionless variable 7' = (ImigWc)"^ intro- 



duced in Sec. |IIE| . The typical value of 7' is of the order 
of 7. We want to estimate the probability density ^(7') 
of large 7', 7' » 7. According to Corollary ^, large 
7' may occur only if Ok ^ Ak for some fc's. Suppose 
for a moment that the number n, which is the small- 
est of such fc's, is larger than two. Let us estimate how 
large the corresponding a„ should be to produce a given 
value of 7'. It is easy to see that (a„/ao)"'^/" ^ yujc/ao 
is needed. In view of Eq. (B£), this can be written as 
(a„/A„)^/" ^ 7'/7- Oi^ the other hand, the probability 
of having a large a„ is, roughly speaking, exp(— a^/yl^). 
This reasoning shows that it is extremely inefficient to 
produce large 7' by boosting a„ with n > 2. Therefore, 
the asymptotic behavior of function P(j') can be ob- 
tained by assuming that all a„ in (Bq) with n > 2 have 
typical values and therefore, can be omitted, while ao, 
ai , and 02 may be large. (This is the smallest number of 
terms needed for the convergence of the remaining inte- 
gral) . It is also quite straightforward to verify that large 
ai's are not as efficient for producing small Imtg's as the 
large values of ao and 02- Hence, ai may be omitted as 
well. Thus, we arrive at the estimate 



1 



JUJc 



> 



ds 



ao + 025 2^0002 



which is the same as 



> 



I'uic ^/2\vv"\ 

To obtain the asymptotic behavior of the distribution 
function P(7') we can replace the inequality sign by 
equality, and so 



00 00 



P{j')~ I dv I dv"5 (7'- ^/^^ ) Vvoh{vy), 







where Prob (u, v") is the joint distribution function of 
V and v" . Neglecting the terms important only for the 
pre-exponential factors, we arrive at 



P(7') ex Prob 



{■KUJcl') 



l\2-\ 



2v, 



where Vs is the saddle-point defined by the equation 
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dv, 



■InProb 



Vs,~- 



(tTWcT') 



l\2-\ 



2v, 



= 0, 



Equation (Bl) can be obtained from here taking into ac- 
count the Gaussian decay of the joint distribution func- 
tion, InProb (w,w") (fMo)^ - (w"M2)^ for v ■> Aq 



and 



Eq. (Bl 



7 ^. 



•d and 



^ ^2. In this simphficd model Vs 
1) follows. 

Finally, let us calculate the numerical factor A, which 
enters Eq. (Bl), for the case of C{r) given by Eq. (C2). 
First we express v and v" in terms of the derivatives of 
Uq. For the sake of notation simplicity we will drop the 
subscript "0" of Uq. In addition, we denote the partial 
derivatives with respect to px by subscript x and with 
respect to Pyhyy. It is easy to see that 



m + U^y 



and 



mujc 



J J _o Tj _o_ 



m 



IJl 



v{p). 



Consider the list 



-\U,Ux,Uy,Ux 



< Uyyy} of U and 



its derivatives arranged in lexicographical order. We will 
refer to the members of this list as uq,ui, . . . ,ug, respec- 
tively. After some algebra, one can derive the following 
results, 



and 



v"v^ 



— —u\ u\ u\ -I- ^2 Wg u\ + u\ui-\- U2 Uq Ui 
-|- 2 U2 Us Wi U4 — 6 U2 "5 U^ W4 -|- 8 1*2 Ui U4 

— 2 U2 ui wg — M2 uf Us + u^us — u\ U5 M3 
-t- 2 U2 U5 u^ U3 — U5 U]^ M3 — 6 u\ ui U4 u^ 

+ 2u2Ui U4 U3 + U2 "^^3 ~ U2 Ui U3 -f U2 U7 



u\ u\ U-! 

ululuQ. 



2 U2 Ui U7 -|- M2 Ui Uq 



Using the fact that Ufe's are Gaussian random variables, 
we obtain 



{u,K~^i 




(B7) 



where iiT is 10 x 10 matrix with elements 

This matrix can be expressed in terms of derivatives of 
Cq{p) at the point p = 0. 



The application of the saddle-point method to the in- 
tegral in Eq. (B7) leads to the estimate 



P(7')oce-™"^i("\ 



where 



Pi(m) = --(«, X-1m), 
and the minimum is sought under the condition 



P2{U) = ^\V^\ 



7r7 ujc 



In fact, we exercise another refinement requiring that 
uo = 0, which reflects the fact that we are interested 
in contours near the percolation level. 

The formulated minimization problem can by solved by 
Lagrange's multiplier method, which requires the mini- 
mization of the form Pi — XP2. It is easy to see, however, 
that it can in its turn be reduced to the (unrestricted) 
minimization of the function Pi{u) — P2{u). Let u* be 
the solution of the latter problem, then the number A 
appearing in Eq. (Bl) is given by 



2 P|(«,)- 



Our result for C from Eq. (CI) is 
A = 4.998. 



(B8) 



APPENDIX C: REALISTIC RANDOM 
POTENTIAL 

It has been suggested that a good model for the ran- 
dom potential really existing in GaAs devices is the fol- 
lowing one: 



C{q) = SnW^d^e 



-2qd 



or equivalently, 

C(r) 



M/2 



(l+r2/4d2)3/2 



(CI) 



(C2) 



Equations (CI) and ( |C2| ) corresponds to the potential 
created in the plane of the two-dimensional electron gas 
by randomly positioned ionized donors set back by an 
undoped layer of width d. The amplitude of the poten- 
tial has the following relation to the parameters of the 
heterostructure 



^2 ^ nn.ie^aBY 



d2 



(C3) 
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where iii is the density of the donors and qb is the ef- 
fective Bohr radius of the electron gas. Equation (|ci^ ) 
applies provided d ^ ub- The random potential can be 
considered a Gau ssian random potential if riiCp ^ 1. 

Using Eq. ( |ClD for the bare potential, one can also 
obtain the real-space correlator Co{p) of the averaged 
potential. Let Re ^ d, then the following relations hold, 



Co{p)=^W^(l-^], 0<p<d, 



where 



W' 



AR^d 



W 



p^ARl - p2 



, p^ d and 2Rc — p ^ d, 






^{Ry,t) = f^ CMmURy)e-<'''-^>^\ (D2) 



M=0 



where function $5j(i?y), given by 



*m(^; 



y- /„»..,,. ... ^'^^m[^ 



VWM!T(7r)i/4 



represents the unperturbed wavefunction of M-th LL 
(Ti-M is the Hermite polynomial). Using this expres- 
sion, one can find the matrix elements UM,M+k{t) = 
{M\U{t)\M + k). If M is large and |fc| < M, then it 
is sufficient to use the quasiclassical approximation [cf. 
Ref. H, Sec. 51] 

UM,M+k{t) - Uk[po{t)], 



Wo = W 



2d 
nRc' 



Note that the 1/p decay of Co{p) for d ^ p <C -Re is a 
universal feature of Co (p) ■ 



APPENDIX D: CALCULATION OF 
QUASICLASSICAL TRANSITION AMPLITUDES 



To derive Eq. (3.3) we start with quantizing the clas- 
sical Hamiltonian (2.2). The result is 



where Uk is the Fourier coefficient defined by Eq. (p.q). 
With the help of this approximation, the equation for the 
expansion coefficient cm can be written as follows 



«?i— 1— = ^ CM+kUk[po{t)]e~ 

k— — oo 



■ikcJct 



It has the solution 

27r 



CM 



de 

2^ 



\o{0)e 



-iMB 



exp|^5fe(t)e^^^M 



H = 



-{Rl + i?,2) + U{p^ + R^,Py -f Ry), 



where Ao(0) depends on the initial conditions at i = io 
and 



where hats are used to denote the operators, viz.. 

ORy apy 



with / = yj%/mujc being the magnetic length. However, 
since the guiding center motion is slow and quasiclassi- 
cal, we can treat px as a classical dynamic variable with 
the equation of motion ( 2_^ ) and similarly for py . 

As in Sec. |l[ we replace the exact trajectory p{t) by 
the "unperturbed" one, po(t). Everything, which was 
said in Sec. Q about the validity of such an approxima- 
tion, applies here as well. 

The Schroedinger equation 



in-'^{Ry,t) = H^{Ry,t), 



(Dl) 



where from now on 



H 



{Rl + Rl) + U[poAt) + Rx,Pov{t) + Ry 



describes the evolution of the cyclotron motion under the 
influence of ti me- dependent perturbation U{t). The so- 
lution of Eq. (Dl) will be sought in the form 



Skit) 



Uk{t')e 



'^"-''"^"^'dt'. 



to 



To elucidate the structure of this solution note that 



formula (D2) can be rewritten in the form 

^{Ry, t)^J2 ^» W *«(^2'' ^) e"^"""*, 



where the relation of $„(i?y,i)'s to $°j(i?y)'s is as fol- 
lows 



M •^ 



e'(n-M)eg^p 



y^ Uk(f)_^~ik{i^at+e) 



fc#0 



khujr 



The new expansion coefficients 6„'s are given by 

27r 



&n = y^A(0)e-"^exp 





E 



Afc(to,t) ^_^feg 
khujc 



(D3) 



where initial conditions now enter through function X{6) 
and Ak{to,t) denotes the following integral 
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Afe(to,0= Uk[po{t')]e 



-ikujct' 



dt'. 



Functions ^n{Ry,t) represent the "instantaneous" LL 
functions at a given point po{t) on the drift trajectory. 
They are the eigenf unctions of H with a "frozen" value 
of po- The corresponding eigenvalues En, however, turn 
out to be time-independent En — ("- + 5) ^^^c + Uq- The 
direct physical meaning have the transitions between the 
instantaneous states $„ not between the unperturbed 
states ^%[- It is the former transition amplitudes we 
are going to calculate (see a similar discussion in Ref. HI, 
Sec. 41). 



After Eq. (D3) is obtained, we can choose any ini- 
tial conditions, for instance, X{6) — e*^^ such that 
bn{to) — Sn,N- In this case bN+kit) gives the desired 
N — > N + k inter-LL transition amplitude, i.e., A^ m+u 
[Eq. (H]. 
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